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Abstract 

The dominant contribution to H~tb production at the LHC is the gluon-gluon 
fusion parton subprocess. We perform for the case of the complex MSSM a complete 
calculation of the NLO electroweak contributions to this channel. The other small 
contributions with quarks or photon in the initial state are calculated at tree level. 
The results are improved by using the effective bottom-Higgs couplings to resum the 
leading radiative corrections. We find that, beyond these leading corrections, the 
NLO electroweak contributions can be still be significant. The effect of the complex 
phases of the soft-breaking parameters is found to be sizeable. 



1 Introduction 



Charged Higgs boson production in association with a top quark is the dominant mech- 
anism in charged-Higgs searches at the LHC The leading order (LO) tree-level diagrams 
involve a gluon and a bottom quark in the initial state. The calculation of the cross section 
can be performed in two ways, by using four- or five-flavor schemes. In the 4-flavor scheme 
(4FS), the bottom density is zero and the leading contribution is gg — > H T tb whose total 
cross section contains large logarithm ~ In fip/rrib, where the factorization scale fiF is of 
the order of the charged Higgs mass. This correction arises from the splitting of a gluon 
into a collinear bb pair. In 5-flavor scheme (5FS) the bottom density is non-zero and the 
leading contribution is gb — > H T t. The large collinear corrections are resummed to all 
orders and are included in the bottom distribution functions. The two schemes should 
give the same result for the total cross section if the calculations are done to a sufficiently 
high order in perturbation theory. A comparison at next-to-leading order (NLO) has been 
done in [TJ. The results of the two schemes are consistent within the scale uncertainties, 
with the central predictions in the 5FS being larger than those of the 4FS [TJ. 

From an experimental point of view, the two final states H T t and H T tb can be sep- 
arated by requiring b tagging. For a heavy charged Higgs boson (M^ > m t ) decaying 
into tb, the signal contains 3te for the former and 46s for the latter. In general, the ad- 
dition of a bottom quark to the final state reduces the signal rate, but the background 
is also lowered. The study in [2] (see also j3] and references therein) shows that a good 
signal-to-background ratio can be achieved by imposing 4 6-tags and suitable cuts if M H ± 
is significantly larger than m t . This study, however, is based on LO predictions and the 
large tan/3 (the ratio of the two vacuum expectation values of the two Higgs doublets) 
enhanced corrections to the bottom-Higgs couplings are not taken into account. Those 
large corrections, which can be resummed and easily included to the LO results by using 
the effective bottom-Higgs couplings, can significantly change the signal cross section, in 
particular for larger values of tan f3. It is therefore important to know the quality of this 
approximation and to have some idea about the remaining higher-order uncertainty. A 
comparison with the full NLO results is needed. 

In the Minimal Supersymmetric Standard Model (MSSM), NLO corrections to charged 
Higgs production in association with heavy quarks at the LHC have been studied to some 
extent. For H T t production, both QCD and electroweak (EW) NLO corrections have been 
calculated (IHE], and some higher-order QCD corrections in 0(10] ■ For exclusive H T tb 
production, the QCD corrections have been considered in [TJ[TT] . The EW corrections are 
missing. All those studies assume that the soft-breaking parameters are real. 

The purpose of this paper is to provide and study the EW corrections to the exclusive 
H~tb production at the LHC for heavy H ± (with > m t ). The tagged bottom quark 
is required to satisfy the kinematic constraint: 

p T:b > 20 GeV, \ Vb \ < 2.5, (1) 

where px,b is the transverse momentum and 775 is the pseudorapidity. The cross section 
after cuts is still considerable. Our study is done in the MSSM with complex parameters 

lr The computer code can be obtained from the authors upon request. 
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(complex MSSM, or cMSSM). The impact of important phases on the cross section will 
be quantified. It turns out that this effect is not small. 

The paper is organized as follows. Section [2]is devoted to the tree- level study, including 
the issue of effective bottom-Higgs couplings. The calculation of the NLO EW corrections 
to the process gg — > H~tb is done in Section[3j Numerical results are presented in SectionH] 
and conclusions in Section [5j 



2 Leading order consideration 

At tree level, the gg contributions of order 0(a 2 a) are dominant. Other contributions 
of the same order arising from qq (q is a light quark) annihilations are much smaller, 
since they involve only the s channel diagrams which are suppressed at high energy and 
the quark density is smaller than the gluon one at the LHC We will, however, include 
those contributions at tree level. It is noted that the qq annihilations give also O (a 3 ) 
contributions coming from tree-level EW Feynman diagrams. These small channels are 
neglected in our calculation. 

We assume the 5FS with b tagging. The three classes of subprocesses of order 0(a 2 a) 

are 

g + g -> H + t + b, (2) 
q + q ->■ H +t + b, (3) 
b + b H + t + b, (4) 

where q = u,c,d,s. The first two channels have been calculated in p^[^ rTTHT3] . The 
corresponding Feynman diagrams of those subprocesses are shown in Fig. [TJ The last 
process is expected to be small and will be shown to be numerically irrelevant. It should 
be noted that the bb annihilation containing the collinear splitting b — > bg* is suppressed 
by the px,b cut. 

There exists also a contribution of order 0(a s a 2 ) arising from the photon-induced 
process, 

g + 7 ->■ H + t + b, (5) 

according to the Feynman diagrams depicted in Fig. [2J Compared to the gg fusion, a new 
EW splitting 7 — > H + H~ appears. This splitting leads to contributions increasing with 
decreasing M^. Although the $7 cross section is larger than the one from bb annihilation, 
it turns out to be negligible as well, as we will show in our numerical analysis. The small 
77 fusion contribution of 0(a 3 ) is neglected. 

All tree-level diagrams involve the Yukawa couplings of the charged Higgs bosons to 
the top and bottom quarks, which read as follows, 

ie ( m t 
X b tH+ = -7= — — 7 a P L + m b t&n (3P R 



y/2s w M w \tan/3 



2 



(a) 




Figure 1: The QCD tree- level diagrams: (a) for the gg — > H tb subprocess, (b) for the 
qq H~tb subprocesses (q = u, c, d, s) and (c) for the bb — > H~tb subprocess. 



where Pl,r — (1 T 7s)/2, sw = sin^^y- It is known that these couplings can get large 
Standard Model (SM) QCD, SUSY-QCD andEW corrections. The SM-QCD corrections 
are absorbed by the replacement nib — > w^° R (/ii?) with fiR being the renormalization 
scale, i.e. the running quark mass is used. The universal SUSY-QCD and EW corrections 
are resummed via the quantity The exact definition of m® R (fiii) and A& are given 
in [H]. We just want to emphasize here that the quantity A;, is proportional to tan/3 
and depends on the mass of the SUSY particles. Including these corrections, the effective 
bottom-top-Higgs couplings read 



V2s w M w V.tan/3 / 

\lH- = fmptan/3AgP L + -^-P fi y (7) 

V2s w M w V tan (3 J 



where 



A 



1 - A 6 /(tan/3) 2 
'' 1 + A 6 



The top-quark mass is considered as the pole mass which is an input parameter in our 
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(*) (*) (*) 

Figure 2: The tree- level diagrams for the #7 — > H~tb subprocess. The star means that 
the graphs with the two incoming particles interchanged also contribute. 



calculation. In the explicit one-loop calculations, we will have to subtract the EW part of 
the Ab correction has already been included in the (IBA) tree-level contribution, to avoid 
double counting. This can formally be done by adding the following counterterms 



H+ 



5m 



5mf" 



m 



DR 



1 + 



m 



DR 



1 



(tan/3) 2 
1 

(tan/3) 2 



(Am 



SEWn* p 



Aml EW P L 



(9) 



to Srrib in the corresponding bottom-Higgs-coupling counterterms, as listed in Appendix B 
of [13]. The definition of Amf w is also given in [ T4 j. 

To quantify the effect we define the improved Born approximation (IBA) where 
the effective couplings in Eq. (J7|) are used. The LO cross section is computed with the 



tree-level couplings in Eq. §6§ with m;, 



m 



DR 



At the end, from the various partonic cross sections, either at LO or IBA, o^o/ibA' we 
obtain the corresponding LO and IBA hadronic cross sections in the following way, 



LO/IBA 



T — 

1,3 



dxidx 2 [F[(x 1 ,ii F )Ff(x 2 ,n F ) ^lo/iba 1 
+ (1^2)], 



a s a,a s a , fx R ) 



(10) 



where = (q,q), (b,b), (g,g), (p, 7); F[(x,^f) denotes the parton distribution function 
(PDF) of parton i at momentum fraction x and factorization scale fip. We use the 
MRST2004qed set of PDFs |17], which includes the photon density in the proton. 



3 NLO electroweak contributions to gg — > H tb 

In this section we discuss the NLO EW contributions to the gg — > H~tb subprocess. 
These corrections are of order 0(a 2 s a 2 ). Other corrections of the same order arising from 
the remaining subprocesses in Eq. (|3J), Eq. (J3J) and Eq. (JS]) are much smaller and will be 
neglected. 

The NLO EW contributions are composed of a virtual part and a real part. The virtual 
part comprises the contributions of bottom-quark and top-quark self-energies, of trian- 
gle, box and pentagon diagrams, and of wave-funtion corrections. For illustration, some 
generic classes of self-energy and vertex diagrams including the corresponding counter- 
terms are shown in Fig. [31 The box and pentagon diagrams are UV finite; a representative 
sample is depicted in Fig. HI 
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Figure 3: Representative one-loop self-energy and vertex diagrams. The shaded regions 
are the one-particle irreducible two- and three-point vertices including the counterterms. 
G denotes the W Goldstone bosons. 




Figure 4: Representative diagrams with irreducible four- and five-point vertices. S denotes 
a Higgs or Goldstone boson, V an electroweak gauge boson, q a squark, and x a chargino 
or neutralino. 
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The virtual part contains UV divergences, soft, and collinear singularities. The UV di- 
vergences are cancelled by renormalization, which requires the choice of a renormalization 
scheme. We follow the renormalization procedure, a hybrid of on-shell and DR schemes, 
described in [TH[T6]|18]. The calculation is done by using the technique of constrained 
differential renormalization (CDR) [19] which is, at one-loop level, equivalent to regu- 
larization by dimensional reduction [20l[2TJ . The on-shell scheme is used for the fermion 
sector, the fine-structure constant, and the charged Higgs-boson mass. The charged Higgs 
field and tan /3 are renormalized in the DR scheme. Since the charged Higgs boson appears 
in an external state, wave-function renormalization is required including also the mixing 
of H~ with W and charged Goldstone bosons (Fig. [3]). 

To make the EW corrections independent of mm./ from the light fermions / 7^ t, we 
use the fine-structure constant at M z , a = a(M z ) as an input parameter. This means 
that we have to modify the counterterm according to 

5Z a(M z ) = §z a(0) _ iAQ^Mf), 



Aa(M| 



<9E 



AA 
T 



dk 2 



Re (M|) 



Ml 



with the photon self-energy from the light fermions only to avoid double counting. In the 
calculation of EW corrections, the couplings in Eq. fl6]) are used. In order to obtain UV 
finiteness, we have to set, as done in [13], ?n& = m^^fl) everywhere. 

We classify the virtual part into two gauge-invariant groups. The first group consists 
of one-loop diagrams contributing to the process 

g + g ->■ H + H + * ->■ H~ +t + b, (12) 

where the virtual H + * can be on-shell, see Fig. [3] and Fig. |4] (box diagrams). The second 
group is the remainder, which is free of resonating propagators. The first group is UV and 
infrared finite since the channel g+g — > H~+H + * does not occur at tree level. Because the 
intermediate H + * can be on-shell, special care has to be taken for the numerical integration 
over the phase space. The resonance propagator reads (zero-width approximation) 

^ = ? - 4 + u = pv (^W) - M[f - Ml±) ' (i3) 

where PV denotes the Cauchy principal value. The principal- value part can be calculated 
by imposing a small cut on q 2 around the pole. The contribution from the S function part 
is nonvanishing because the imaginary part of the on-shell propagator can multiply by 
the imaginary part of the loop integrals, hence the corresponding one-loop amplitude can 
interfere with the tree-level amplitude. We have checked that this contribution is indeed 
nonzero, but small. A naive calculation taking into account only the principal value part 
would lead to an incorrect result. For practical purposes, a better method is introducing 
a small width in the resonance propagator, 

AH± = q 2 -M 2 H± l +l M H ,T H± - (14) 
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We have checked that the result is practically independent of the small values of the 
width and agrees with the sum of the principal value and 5 function contributions. We 
also notice that this method gives smaller integration error. As will be shown in the 
numerical study, the effect of the H~H + * production mechanism is small at the cross 
section level, but is of importance for differential cross sections. 

The real EW corrections arise from the photonic bremsstrahlung process, 

g + g -> H +t + 6 + 7, (15) 

with the corresponding Feynman diagrams shown in Fig. [51 This contribution is divergent 
in the soft limit (p° — > 0) and contains quasi-collinear corrections [25J proportional to 
a\og(ml/E^), E b being the 6-quark energy, in the limit pbP^ — > 0(ml). The 6-quark mass 
is used for regularization and to separate the singular terms. A fictitious photon mass (A 7 ) 
is used for regularization of the soft singularities. If we consider the total cross section, 
i.e. without applying the cuts in Eq. fll}), the soft and quasi-collinear singularities cancel 
completely in the sum of the virtual and the real contributions, according to the Kinoshita- 
Lee-Nauenberg theorem [221123] • This requires that we have to use = mP R (/i#) as in 
the virtual amplitudes. If the cuts in Eq. ([T]) are imposed then the soft singularities 
still cancel, but the quasi-collinear singularities do not, since the cuts requiring bottom- 
photon separation are not collinear safe. In this case, some quasi-collinear singularities 
remain and are regularized by the bottom mass. Those left-over singularities can be 
separated, as discussed below If a sufficiently collinear 6-photon system is recombined 
before applying cuts then the quasi-collinear singularities cancel, but the result will depend 
on the recombination parameter. As done in the previous study for the NLO QCD 
corrections pQ, we assume in this paper bottom-photon separation, and hence no photon 
recombination is applied. 

The dipole subtraction method [24T426] is used to extract the singularities from the real 
corrections and combine them with the virtual contribution. The subtraction method for 
doing the phase-space integration for the radiation process Eq. ( |15p arranges the integral 
in the following way, 

cw = / [da Icsd 6(p b ) - da suh 6(p b )} + / da suh 6(p b ). (16) 

J A J A 

The subscript 4 refers to the 4-body final state including the radiated photon, 6 is a func- 
tion to impose the kinematical cuts defined in Eq. (Q]), p& is a function of pi with % = H~, 
t, b, 7, with the definition given in [24TT26] . The subtraction function da sn h0(pb) has to be 
chosen such that the first integral is finite and the second one can be partially analytically 
integrated over the singular variables. The function <i<7 su b has the same singular structure 
as da Tea _\ pointwise in the phase space. There are two ways to deal with the cut function, 
i) We require that (the pseudorapidity cut is neglected to simplify the discussion) 

0{pb) -> 0(p b ) or p b -)■ p b (17) 

in the singular limits (the soft limit is trivially satisfied), which implies that 6(pb) is 
not collinear safe, so that the first integral is soft and (quasi-) collinear finite. All soft 
and quasi-collinear singularities are contained in the second integral. All soft and some 
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Figure 6: Partonic cross sections as functions of the center-of-mass energy v£ in the left 
panel and of the cutoff parameter 6 S in the right panel are presented for the gg — > H~tb r y 
process. The results are obtained by using the phase space slicing method (PSS) and the 
dipole subtraction method (SUB, straight line) in the CPX scenario. The photon mass 
regulator In A 7 is set to zero. 



quasi-collinear singularities are canceled in the sum with the virtual contribution. The 
leftover quasi-collinear singularities, regularized by nib, can be factorized and separated. 
A detailed procedure including the definition of pb is described in [22] . A consequence of 
the condition (TTTj) is that, in the calculation of the first integral, we can set nib = in the 
kinematics (but not in the Yukawa couplings). 

ii) We require that the cut function 9(pb) is infrared safe as in [2^1 [25] so that the 
sum of the second integral and the virtual contribution is independent of soft and quasi- 
collinear singularities. Specifically, it means that the condition ( 1T7|) is satisfied for the 
soft limit but not for the collinear limit. The first integral, therefore, contains the leftover 
quasi-collinear singularities. Since the result is finite one can do it numerically. In this 
approach, one has to keep nib everywhere. We have implemented both approaches and 
found good agreement for cross section and distributions. Moreover, the result of the 
dipole subtraction method is compared with the one of phase-space slicing method, as 
illustrated in Fig. El In the numerical analysis, we will present the results of the dipole 
subtraction method because the integration errors are smaller. 

The above treatment of kinematical cuts in the dipole subtraction method is also 
applied for the bottom-quark histograms displayed in Section 14.41 

Finally, the hadronic cross section at NLO is written in the following way as the sum 
of the improved Born approximation and the genuine loop and radiation terms for the gg 
subprocess, 

pp _ pp _j_ APP/sQ 

NLO — IBA + U EW ' 1 15 ) 
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with 



Aew 9 = J dxidx 2 F*( Xl , fi F )FP(x 2 , fi F ) A™ w (a 2 s a 2 , p R ) . (19) 

Thereby, is the sum of the virtual and real corrections at the partonic level, as 

discussed above. The IBA part of the cross section results from the various sources at the 
partonic level, 

rrVV _ , -pp/qq j_ vp/bl pp/gi / on \ 

°IBA — IBA "+" °IBA + IBA + °IBA > l zu J 

as discussed in Section [2j 



4 Numerical studies 
4.1 Input parameters 

We use the same set of input parameters as in [H] for the sake of comparison. For the 
SM sector: 

a s {M z ) = 0.1197, a(M z ) = 1/128.926, 

M w = 80.398 GeV, M z = 91.1876 GeV, (21) 

m t = 173.1 GeV, m b (m b ) = 4.2 GeV. 

We take here a s = a^ s (fj,R) at three-loop order [27] . m^mj) is the QCD-MS 6-quark 
mass, while the top-quark mass is understood as the pole mass. The CKM matrix is 
approximated to be unity. 

For the soft SUSY-breaking parameters, the adapted CP-violating benchmark scenario 
(CPX) [281129] is used, 

\fi\ = 2 TeV, |M 2 | = 200 GeV, |M 3 | = 1 TeV, \A t \ = \A b \ = \A T \ = 900 GeV, 
Mq = M b = M = M- L = M E = Msusy = 500 GeV, \M X \ = 5/3 tan 2 9 W \M 2 \. ( " ] 

We set Af = for / = e, /1, u, d, c, s since the Yukawa couplings of the first two fermion 
generations proportional to the small fermion masses are neglected in our calculations. 
With the convention that M 2 is real, the complex phases of the trilinear couplings A t , A b , 
A T and the gaugino-mass parameters Mi and M3 are chosen as default according to 

7T 

4> t = (f> b = (fi T = 01 = 03 = ^ ' ( 23 ) 

unless specified otherwise. The phase of /1 is chosen to be zero in order to be consistent 
with the experimental data of the electric dipole moment. We will study the dependence 
of our results on tan 0, Mh± and 0t in the numerical analysis. 

The scale of a s in the SUSY-QCD resummation of the effective bottom-Higgs couplings 
in Eq. (9) of [2] is set to be Q = {ra bl + m~ b2 + m§)/3. If not otherwise specified, we 
set the renormalization scale equal to the factorization scale, = /xp, in all numerical 
results. Our default choice for the factorization scale is //po — ( m t + M H ±)/3. 
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Our study is done for the LHC at 7 TeV, 8 TeV and 14 TeV center-of-mass energy 
(a/s). In the following we show the dependence of the cross section on tan/3, M H ± and 4>u 
and various differential distributions for the default parameter point. Since the results of 
the different center-of-mass energies look quite similar in shape and differ mainly by the 
magnitude of the cross section, our discussion essentially applies to all displayed cases of 
the total energy 

4.2 Calculations and checks 

The results in this paper have been obtained by two independent calculations. We have 
produced, with the help of FeynArts-3.4 [301I3T] and FormCalc-6 . [21], two different 
Fortran 77 codes. The loop integrals contain five-point tensor integrals up to rank three, 
and four-point tensor integrals up to rank three. The pentagon integrals are reduced to 
the box integrals by using the reduction methods in [32TI34"] . The two-, three- and four- 
point tensor integrals are further reduced to the scalar integrals by using the Passarino- 
Veltman reduction method [35]. The loop integrals are evaluated with two independent 
libraries, LoopTools/FF [3"EH3"8"] using the five-point reduction method of [32J and our 
in-house library Looplnts using the five-point reduction method of [51] ■ The latter uses 
the method of [391 - HT] . treats all the internal masses as complex parameters and has an 
option to use quadruple precision, on the fly, when numerical instabilities are detected. 
The phase-space integration is done by using the Monte Carlo integrators BASES [12] and 
VEGAS [13] . The results of the two codes are in good agreement. On top, we have also 
performed a number of other checks: 

For the process gg — > H~tb, we have verified that the results are QCD gauge invariant 
at LO, IBA and NLO. This nontrivial check, which can detect a bug in the Feynman rules 
and in the tensor reduction procedure, can be easily done in practice by changing the 
numerical value of the gluon polarization vector e^(p, q), where p is the gluon momentum 
and q is an arbitrary reference vector. QCD gauge invariance means that the squared 
amplitudes are independent of q. More details can be found in [11]. The common checks 
of UV and infrared fmiteness are done for the NLO calculations. 

4.3 LO, IBA and NLO cross sections 

The LO, IBA and NLO cross sections as functions of tan/3, the phase </> 4 , and M H ± are 
shown at 8 TeV and 14 TeV on the left and right columns of Fig. [71 respectively, and at 
7 TeV on the left column of Fig. [8] The relative corrections 5, with respect to the LO 
cross section, are defined as 5 = (ct IB a/nlo — olcO/olo 

We first study the effects of A;, resummation in the effective bottom-Higgs couplings. 
For small values of tan/3 the left-chirality contribution proportional to m t / tan/3 is dom- 
inant while the right- chirality contribution proportional to mb tan /3 dominates at large 
tan/3. The cross section has a minimum around tan/3 = 8. The effect of resumma- 
tion is best understood by comparing the phase-dependence plot and the others. The 
important point is that Aj, is a complex number and only its real part can interfere with 
the LO amplitude. Thus, the Af, effect is minimum at <p t = ±7r/2 where the dominant 
contributions are purely imaginary and is largest at <p t = 0, ±7r. The phase- dependence 
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plot shows that the Amjp w effect can be more than 10%. From the tan/3 dependence 
plots where A b is mostly imaginary we see the effect of order C(A^), which is about —15% 
at tan /3 = 10. 

We now turn to the NLO cross sections, which include the complete EW corrections 
to the process gg — > H~tb. Fig. [7] also contains the effect of the H~H + * resonance 
mechanism, which is almost invisible at the cross section level. The NLO cross section 
depends strongly on <p t . The IBA results are always closer to the NLO values rather than 
to the LO ones. In particular, for the phase dependence, the IBA shows the qualitative 
features of the NLO prediction while the LO cross section is a constant. After subtracting 
the Am^ EW corrections, the remaining NLO EW contributions are still sizeable. The 
relative correction increases with tan/3 and M^± for the default value <pt = 7r/2; for fixed 
default values of tan/3 and M#±, it is maximum (about 40%) at <p t = 0. As an aside, we 
remark that the IBA and NLO EW effects for the process gg — > H~tb are similar to the 
ones found in the bb — > W~H + study [T4] . 

Table 1: The total cross section in fb for pp — > H tb including the IBA of the four subprocesses and 
EW NLO corrections to gg — > H~tb at y 7 ! = 14 TeV. The charged Higgs-boson masses are given in GeV. 
The numbers in brackets show the integration uncertainty in the last digit when they are significant. 
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Table 2: Similar to Table [TJ but for y/s = 8 TeV. 
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600 


0.18947(5) 


0.00444 


0.00056 


0.00315 


-0.02328(8) 


0.17435(10) 



Table [T] shows separately the IBA results for the various subprocesses together with 
the NLO EW corrections to gg — > H~tb at a/s = 14 TeV for different values of M#± and 
tan/3. Similar results are presented in Table |2] and Table [21 but now for y/s = 8 TeV 
and y/s = 7 TeV, respectively. We observe that the gg contributions are dominant; they 
contribute more than 90% (83%) of the total IBA for ^Ts = 14 TeV (y/s = 7 TeV). The 
contribution of the bb channel is below 1%; the #7 channel contribution is slightly larger. 



11 



Table 3: Similar to Table Q] but for V* = 7 TeV. 



tan/? 


M H ± 


„pp/gg 

°IBA 


jpp/qq. 

"IBA 


pp/bb 
°IBA 


pp/g-y 
°IBA 


App/gg 


all 


5 


200 


5.3652(9) 


0.9885 


0.0311 


0.1058 


-0.2049(6) 


6.286(1) 


10 


200 


3.5138(6) 


0.6551 


0.0205 


0.0695 


-0.3552(5) 


3.9037(8) 


20 


200 


6.085(1) 


1.095 


0.035 


0.119 


-1.367(2) 


5.967(2) 


10 


300 


1.2570(2) 


0.0974 


0.0058 


0.0224 


-0.1292(2) 


1.2534(3) 


10 


400 


0.5164(1) 


0.0242 


0.0019 


0.0089 


-0.0544(1) 


0.4971(1) 


10 


600 


0.10583(2) 


0.00268 


0.00027 


0.00191 


-0.01295(2) 


0.09774(3) 



The NLO EW contributions are comparable in size to the qq contributions, but with the 
opposite sign. 

4.4 Differential distributions 

We now consider the differential distributions of various kinematical variables, in the IBA 
and including the NLO EW corrections. The relative correction is defined with respect 
to the IBA differential cross section, 5 = (da^LO — dam a) / dam a- All results are shown in 
the right column of Fig. [8] and Figs. (0 [ID], ITT]) . 

The effect of the H~H + * production mechanism is best seen in the right column of 
Fig. [HI The tb invariance mass distribution shows the singular pole structure at M t i = 
Mh± = 200 GeV if this channel is included. This effect is also visible in other distributions. 

Distributions for the individual particles separately are shown in Fig. [91 for y/s = 
14 TeV, and in Figs. f Til3]rrT]) for the lower energies 8 and 7 TeV. The results are very 
similar and differ essentially in the absolute size of the cross section at the different 
energies. 

For the charged Higgs boson, the relative correction is negative, decreases with pt(H~) 
and has a minimum (about —10%) at the central rapidity. 

For the top quark, the behavior of the px distribution is similar to the one of the 
charged Higgs boson. The EW corrections are negative and decrease with p?, consistent 
with Sudakov corrections a log n (p T /M^ v ) with n = 1, 2. For the rapidity distribution, the 
relative correction is rather flat (about —10%) in the region \y t \ < 3. 

The distributions of the bottom quark are quite different from the ones of the heavy 
particles. At tree level (see the IBA curve), the cross section is larger at low pt due to 
collinear bottom-quark radiation off gluons. The relative correction increases and reaches 
the maximal value at pt ~ 0.3 TeV and then follows the trend of decreasing with p? 
as for the other particles. This behavior can be explained by the interplay between the 
leading weak Sudakov correction a\og 2 {p T / M^) and the QED quasi-collinear correction 
a \og{m1/p\) from photon radiation off the bottom quark. The latter is more important at 
low pt while the former dominates in the high energy regime. For the rapidity distribution, 
the correction is smallest in the central region. 
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5 Conclusions 



In this paper we have studied the production of charged Higgs bosons in association with 
a top quark and a tagged bottom quark at the LHC in the context of the complex MSSM. 
Cuts on the transverse momentum and rapidity of the bottom quark are applied. At tree 
level, the gg fusion is dominant among various subprocesses with quarks or photon in the 
initial state; for this parton process, the NLO EW corrections have been calculated and 
discussed. 

Since the tree-level amplitudes are proportional to the top-bottom-Higgs coupling, 
we have examined the effective-coupling approximation and compared it to the full NLO 
result. The dependence of the cross section on tan j3, Mh± and the phase (f> t of the trilinear 
coupling A t has also been studied. 

Numerical results have been presented for the CPX scenario. The production cross 
section shows a strong dependence on tan/3, M H ± and 4>t- Large production rates occur 
for small tan /3, small Mjj± and phases <f>t around ±7r. At LO, the cross section increases 
strongly with large tan f3. This behavior is, however, significantly reduced when NLO 
corrections are included. An interesting feature is the <\> t dependence: while the LO cross 
section is just a constant, the IBA and NLO results show a strong dependence with a 
minimum at <pt — 0. 

We have also presented various differential distributions of the final state particles, 
where the NLO EW corrections are usually negative. 
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Figure 7: The cross section for pp — > H~tb as function of tan /3, (j) t , Mh± in various steps 
of approximation LO, IBA, and NLO, for yfs = 14 TeV (left) and y/s = 8 TeV (right). 
The lower part of each panel shows the relative corrections to the LO results. Also shown 
is the NLO result without the H~H + * production mechanism. 
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Figure 8: Left column: the LO, IBA and NLO cross section for pp — > H~tb as function of 
tan /3, 4>t, and Mh±, for yfs = 7 TeV. Right column: the IBA and NLO EW invariant mass 
distributions of the tb system for pp — > H~tb at 14, 8 and 7 TeV. The lower panels show 
the relative corrections. Also shown is the NLO result without the H~H + * production 
mechanism. 
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Figure 9: The IBA and NLO distributions of transverse momentum (left) and rapidity 
(right) for H~, t, b at 14 TeV. The lower panels show the relative corrections. Also shown 
is the NLO result without the H~H + * production mechanism. 
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Figure 10: Similar to Fig. |H]but for y/s = 8 TeV. 
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Figure 11: Similar to Fig. |H]but for y/s = 7 TeV. 
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